Three-dimensional vortex dynamics in Bose-Einstein condensates 
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We simulate in the mean-field limit the effects of rotationally stirring a three-dimensional trapped 
Bose-Einstein condensate with a Gaussian laser beam. A single vortex cycling regime is found for a 
range of trap geometries, and is well described as coherent cycling between the ground and the first 
excited vortex states. The critical angular speed of stirring for vortex formation is quantitatively 
predicted by a simple model. We report preliminary results for the collisions of vortex lines, in 
which sections may be exchanged. 
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The formation and properties of vortices in Bose- 
Einstein condensates are attracting intense current in- 
terest. A number of schemes for producing vortices in 
trapped condensates have been proposed theoretically, 
and the first experimental creation of a topological vor- 
tex in a two-component condensate has been reported jl] , 
in agreement with a numerical model O . Simulations of 
flow of a condensate past a spherical object have shown 
vortex ring solutions |3|] , and numerical studies of vortex 
formation in rotating traps have shown the formation of 
vortex arrays ||]. Very recently the first report has been 
made of the observation of vortices in a single-component 
condensate ||. 

In an earlier paper gj], we showed that a rotationally 
stirred Bose-Einstein condensate, simulated in two di- 
mensions, exhibits a simple single vortex cycling regime, 
a behaviour we described as nonlinear Rabi cycling. In 
the present paper we demonstrate that the single vor- 
tex cycling regime is also found in full three-dimensional 
simulations of condensate stirring (as shown in Fig. [I]), 
in a parameter regime accessible to current experiments. 
The single vortex regime is studied in a variety of trap 
geometries (oblate to prolate) as is the dynamical stabil- 
ity of central vortices. The line character of a vortex in 
three dimensions also gives rise to some qualitatively new 
behaviour that has no analogue in two dimensions. We 
show that in a spherical trap, stirring of the condensate 
can produce two vortex lines which collide and exchange 
sections. 

As in the earlier work, our treatment is based on the 
time-dependent Gross-Pitaevskii (GP) equation, 



[-V 2 + V(r,i) + C|V(r,i)| 2 ] ,/>(r,t), (1) 



in dimensionless units 
V(r, t) is given by 



x 2 +y 2 + X' 



The total external potential 
" z 2 )/4 + W(r,t), where the 



first term represents a cylindrically symmetric trap with 
trap anisotropy parameter A. The second term W(r,t) is 
the contribution of the stirrer, a far-blue detuned Gaus- 
sian laser beam cylindrically symmetric about a line par- 
allel to the z axis, 




FIG. 1. Probability density isosurface (\tj)\ = 10~ ) for a 
condensate stirred as described in the text, with the stirrer lin- 
early withdrawn between t = 7n and t = 8ir. Here t — 94.12, 
and the hole near the centre of the condensate is a vortex 
that has been drawn in from the edge and remains close to 
to the centre of the condensate. Parameters are C = 1000, 
A = Wo = 4, w s — 4, and p s = 2. The angular fre- 
quency of stirring ulf — 0.3 is less than the critical frequency 
LJc = 0.323. 



W(r,t) = W s (t)exp 



\P-Ps(t)\ 

Ws/2 



(2) 



where p is the projection of r into the z — plane, and 
p s is the centre of the Gaussian stirrer in that plane. 

The condensate is stirred in a manner similar to our 
earlier work ||, with the stirrer moving in a circle of 
radius p s at a constant angular velocity ujf. Here, how- 
ever, we begin with the ground state of the harmonic trap 
(with no stirrer present), and minimise transient effects 
by increasing the stirrer amplitude W s linearly from zero 
to a final value of Wq between t — and t = it. With 
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FIG. 2. Phase of ip in the z = plane for the condensate 
shown in Fig. |l|. 



100 



FIG. 3. Condensate angular momentum expectation value 
versus time, for the case of a constantly rotating stirrer (thin 
line), or stirrer linearly withdrawn between t — 7n and t = 8n 
(thick line). Parameters as in Fig. [[]. 



an appropriate choice of parameters, as in Fig. [l|, a single 
vortex with positive circulation is drawn close to the cen- 
tre of the condensate, and for continuous stirring cycles 
in and out. This is exactly the behaviour found previ- 
ously in a two-dimensional system Q], and in fact the 
parameter choice for the single vortex cycling behaviour 
in three dimensions is guided by the same considerations 
outlined there, as we discuss below. Figure ||, where we 
plot the phase of the wavefunction, shows that the cen- 
tral feature of Fig. [j] is a vortex of positive circulation. 
In Fig. H the evolution of the angular momentum expec- 
tation value (which in our dimensionless units is for 
the ground state and 1 for the central vortex state) is 
presented and confirms that for continuous stirring (thin 
line), the vortex cycles in and out of the condensate. 

In contrast to the two-dimensional case, the three- 
dimensional simulation can be unambiguously related to 
a realistic current experimental scenario. The C = 1000 
case shown in Figs. [j]-|3| corresponds to N — 1.8 x 10 4 
atoms of 87 Rb in the \F — 2, nif = 2) hypcrfme state (for 
which we use an s-wave scattering length a = 5.29 nm) 
in a time orbiting potential (TOP) trap with radial trap 
frequency uj r = 2n x 15 Hz and A = w 2 /w r = \/8. For 
this example, our computational units of time and length 
M correspond to 10.6 ms and 1.97 /im respectively, and 
the Thomas-Fermi diameter of the condensate in the ra- 
dial direction of 11.7 corresponds to about 23 /im. The 
stirrer Gaussian 1/e diameter w s = 4 corresponds to a 
Gaussian 1/e 2 beam diameter of approximately 11 /im. 
The required intensity of the stirring laser is detuning de- 
pendent, but can be easily calculated from the light shift 
potential M7 2 /4A, where f2 and A are the Rabi frequency 
and detuning respectively of the atom-laser interaction. 
For example, if the laser beam is 50 nm blue-detuned 
from the 780 nm atomic transition, A = 1.65 x 10 14 s _1 , 
and a laser power of 300 nW is required for the 3 nK 



height stirrer we describe. The first angular momentum 
peak in Fig. [^occurs at 0.26 s, which is much less than the 
lifetime of the condensate. The occurrence of the single 
vortex cycling behaviour in this experimentally accessi- 
ble parameter regime suggests a method of preparing a 
central vortex state; stir the condensate as we have de- 
scribed, then remove the stirrer at the time when the 
angular momentum peak is expected. We illustrate the 
implementation of this method in Fig. [?] with the thick 
curve, for which the stirrer is linearly withdrawn between 
t = 7tt and t = 8tt. Once the stirrer is removed, the vor- 
tex remains in the condensate, circling about and close 
to the centre, and (L z ) is constant. 

The vortex cycling behaviour illustrated by the ex- 
ample in Figs. [l]|] indicates that the two-state model 
described in our earlier work (for the two-dimensional 
case ||]) also has application in three-dimensions. Here 
we briefly outline the modifications required to extend 
the two-state model into three dimensions. We represent 
the condensate by an axisymmetric component <p s , and 
a vortex-like component cf) v e ie which has a circulation 
about the z axis, where <fi s and <p v are real and nonneg- 
ative, so 

ip(r,t) = a s (t)tp s (p,z,n v ) + a v (t)cj) v (p,z,n v )e l6 ', (3) 



with p — \Jx 2 + y 2 . We choose 4> s and <f> v to satisfy the 
following time-independent coupled equations, 



ld_ ( d_ 

p dp \ df 
p 2 + X 2 z 2 



dz 2 



+ C (n s 4> 2 s + 2n v 4> 2 v ) 



(4a) 
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P 2 + X 2 



+ C {n v <p 2 v + 2n s cf> 2 s 



(4b) 



where the eigenvalues (i s and fx v are found by solv- 
ing Eqs. (4b) for a given choice of the vortex fraction 



n v = \a v \ . The quantity n s = \a s \ — 1 — n v . We note 
that when n v = 0, Eqs. (Q) reduce to Eq. ([la]), the time- 
independent GP equation for the ground state, while if 
n v = 1, Eqs. (Q) reduce to Eq. (4b), the time- independent 
GP equation for the first excited vortex. The result- 
ing nonlinear Rabi equations for the two-state system 
[obtained by substituting Eq. (^) into Eq. (|l|)] are then 
identical in form to Eqs. (5-6) of [||, but with the wave- 
functions <p s and <p v now defined over p and z instead of 
r. 

The underlying assumption of our two-state model is 
that the ground state (or more generally, <p s ) is coupled 
primarily to the vortex state of lowest energy with axis 
parallel to the stirrer axis. In two dimensions, there are a 
limited number of excited condensate states in the appro- 
priate energy range that can be coupled to the symmetric 
state by the stirring potential, and we have previously 
shown j(| that the two state model works well within 
a specified validity range. In three-dimensional cylindri- 
cally symmetric traps, many more condensate eigenstates 
exist, and although symmetry considerations limit the 
possible coupling from the symmetric state, a two state 
model places more constraint on the description of the 
system than in the two-dimensional case. We do not aim 
therefore to provide a detailed representation of the stir- 
ring behaviour with our two state model, instead we use it 
primarily to provide conceptual understanding of the cy- 
cling behaviour. Nevertheless we also find it gives a quan- 
titatively accurate prediction of the critical frequency of 
rotation uj c , the stirring frequency which causes a single 
vortex to cycle right to the centre of the condensate. We 
have examined the single vortex cycling behaviour us- 
ing the full GP solution, for a variety of trap geometries 
including the oblate (A = spherical (A = 1), and 

prolate (A = 1/3) cases, with a condensate of C — 1000, 
and a stirring potential of the same size and position as 
in Fig. [l]. We find that the critical rotational frequency 
for single vortex cycling is accurately predicted by the 
two state model as the difference in energy between the 
ground and vortex states of the condensate (modified by 
the finite size of the stirrer, see Eq. (7) of ||). As in 
H , a sufficiently large Rabi frequency is required to dis- 
tort the energy barrier between ground and vortex state 
enough to permit cycling. When the critical frequency 
of rotation is exceeded, additional vortices penetrate the 
condensate. 

It is worth noting that as the z component of the 
trapping potential become weaker (that is, as A becomes 
smaller), then in addition to the vortex state, the stirring 
potential increasingly excites the dipole (centre of mass) 
oscillation of the condensate. This can be understood in 
terms of the Ehrenfest theorem: for our choice of param- 
eters the mean force (—WW) of the stirrer increases as 



the condensate radial size decreases because the conden- 
sate is increasingly concentrated near the steepest part 
of the stirring potential. The centre of mass motion con- 
tributes angular momentum (L z )com to the condensate 
angular momentum, where 



{L z )coM = {x){P v )-{y){P x ). 



(5) 



We find that (L z ) — (L z )com (which is a measure of 
the angular momentum caused by the presence of the 
vortex) for all three trap geometries cycles between 
and 1 when to j is slightly less than lo c (the single vor- 
tex cycling regime) and oscillates about 1 when ujf just 
exceeds u c (and multiple vortices penetrate the conden- 
sate). (L z )com is small for the A = case shown in 
Fig. | 

If the single stirrer used in Fig. |lj is replaced by two stir- 
rers of half the height but on opposite sides of the z axis, 
no vortex cycling occurs, and (L z ) exhibits only small 
oscillations of amplitude less than 0.04. This is easily 
understood from our two state model, because the Rabi 
frequency (see Eq. (6b) of ||) is zero for this stirrer 
configuration, as it is for any stirrer configuration of even 
symmetry. It is worth pointing out that our scenario of 
coherent cycling should be distinguished from vortex pro- 
duction by damping in a trap with a symmetric rotating 
distortion, such as the recently reported experiment [^), 
or numerical damping schemes such as imaginary-time 
propagation 

Our three-dimensional simulations also allow us to test 
(in the mean-field limit) the stability of pure vortex 
eigenstates. Using as an initial state a central vortex 
(an eigenstate of the time-independent GP equation) for 
C = 1000, a perturbation was applied by linearly insert- 
ing and withdrawing (between t — and t = ir) a stirrer 
of the same size as in Fig. [j] but at a fixed position in 
the laboratory frame. We find that for a range of trap 
geometries (A = \/8, 1, or 1/3), the (L z ) = 1 vortex is 
dynamically stable, and remains near the centre of the 
condensate. We can understand this result in terms of 
Rabi cycling. The stationary stirrer (tOf = 0) provides a 
potential which is far-detuned from resonance (ujf = lo c ), 
and thus there is a low probability of transfer out of the 
initial (vortex) state. Indeed if the stirrer remains fixed 
in the condensate after insertion and is not withdrawn, 
we find that the vortex still remains near the centre of 
the condensate and follows a small closed path displaced 
slightly away from the stirrer. By contrast, for these ge- 
ometries, the (L z ) = 2 central vortex eigenstate immedi- 
ately dissociates into two unit vortices when perturbed. 
We note that these stability results refer to dynamical 
stability in the GP limit (that is, T = 0), as opposed to 
thermodynamic stability 

In addition to confirming the basic features of the two 
state model of vortex cycling, our three-dimensional sim- 
ulations also reveal some new phenomena. The vortices 
produced in three dimensions by stirring are line vor- 
tices. In the A = V8 configuration, the lines are nearly 



3 




FIG. 4. Collision of vortex lines and exchange of sections 
in a stirred condensate. View is from a position in the z — 
plane, looking almost down the y — —x line. In (a) t — 6.79, 
(b) t = 6.91, and (c) t = 7.04. Parameters are the same as 
in Fig. ^ except that A = 1, The angular frequency of stirring 
u>f — 0.45 exceeds the critical frequency u> c — 0.399. 

straight and parallel to the stirrer axis, and when two 
lines of the same circulation form, we find they repel 
each other. For a spherical trap (A = 1) however, vor- 
tex lines along any direction are degenerate (in the ab- 
sence of the stirrer) with the vortex eigenstate about the 
z axis. Accordingly, the vortex lines can become more 
curved than in the oblate case, as we find in our simula- 
tions below the critical frequency, where the vortex that 
forms can have appreciable curvature. Above the criti- 
cal frequency, when multiple vortices are drawn into the 
condensate, there is a striking consequence of this poten- 
tial for vortex lines to deform. In a collision between two 
such curved vortex lines, sections of the lines can be ex- 
changed, as we show in the collision sequence in Fig. ^. 
The condensate is stirred somewhat above the critical 
frequency in a spherical trap, and by t — 6.79 two vortex 
lines (both of positive sense {!(]]) have formed. The fig- 
ure shows only the vortex lines (detected numerically), 
which have been shaded light and dark to distinguish 
one from the other. In Fig. ^(a) the dark vortex line ap- 
proaches the light vortex line from behind, and because 
of their mutual repulsion, the light line has developed 
a bulge. In Fig. [|(b) the dark line is now close to the 
light line, causing the bulge to become recurved, with 
two small portions now nearly antiparallel to the dark 
line. In Fig. |](c) the two lines have completed a collision 
and exchanged their central sections. 

In conclusion, we have shown that the single vortex 
cycling regime we discovered previously in simulations 
of two-dimensional condensates is also present in three- 
dimensional condensates, for a range of trap geometries. 
We have related our simulations to current experimen- 
tal configurations, and have demonstrated how the cy- 



cling behaviour could be used to generate a central vor- 
tex state. We have investigated the dynamical stability of 
central vortex states, and have also given initial evidence 
of the rich dynamical behaviour of vortex lines. 
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